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I. INTRODUCTION 

Consider a system of non-interacting passive particles 
that undergo a combination of diffusion and convection 
that is driven by a steady velocity field v = v(r). The 
particle density c = c(r, t) obeys the equation of motion 



dc 
dt 



(1) 



where D the diffusion coefficient. If the velocity field is 
radial, v = v{r) f, and the initial conditions are radially 
symmetric, the density satisfies 



dc dc „ 
at or 



1 dc 
dr 



(2) 



Suppose that the velocity field is inversely proportional 
to the radial distance, v = ^P-. Such a divergenceless 

' 27rr ^ 

field is natural in two dimensions and is generated by a 
point sink (or source) of strength Q in an incompressible 
fluid. This flow fleld has the remarkable property that 
the convection term can be absorbed into the diffusion 
operator by an appropriate shift of the spatial dimension. 
Indeed, the density satisfies 



dc 
di 



D 



d^ 



-Icff 



_! dc 

dr 



where the effective spatial dimension is given by 



Q 

2ttD 



= 2-2A, A = 



A-kD 



(3) 



(4) 



Here A is the Peclet number; the numerical factor (47r)^^ 
is chosen to simplify formulae that follow. 

We can thus interpret convection-diffusion in two di- 
mensions with a radial velocity v oc r~^ as isotropic dif- 
fusion in a space with variable effective dimension. The 
motivation for developing this connection is that there 
are many problems in many-body physics for which the 
spatial dimension is an important determining factor in 
the phenomenology. In diffusion-controlled reactions, for 
example, a large body of work has uncovered the gen- 
eral feature that when the spatial dimension d exceeds a 
critical value dc (that depends on the specifics of the re- 
action), the reaction kinetics has a mean-fleld character, 



a situation where any pair of reactants is equally likely 
to react. Conversely, when the spatial dimension d < dc, 
fluctuation-dominated phenomena arise, such as anoma- 
lously slow kinetics and non-trivial spatial organization 
of reactants, see, e.g., Refs. [E 0, 0, S IE 13 for a review. 
The radial flow system offers a potentially attractive way 
to study the full range of behavior between the disparate 
mean-field and fiuctuation-controUed regimes simply by 
tuning the flow velocity. 

In this work, we therefore exploit the connection be- 
tween convection-diffusion in two dimensions and pure 
diffusion in a space with a tunable effective dimension. 
By this equivalence we can probe diffusion-controlled pro- 
cesses in arbitrary dimension by varying the Peclet num- 
ber. In the next section we emphasize some subtleties 
associated with this mapping. In the following sections, 
we study simple diffusion-controlled reactions in general 
dimensions by exploiting the mapping to radial flow in 
two dimensions. We first consider the influence of an ab- 
sorbing trap on the concentration proflle in Sect. ITTll In 
particular, we analyze the minimal separation between 
the particles and the trap. In Sect. lIVl we consider two 
simple reactive systems, namely, irreversible coalescence 
and irreversible annihilation. 



II. SPREADING OF A RING 

As a warm-up exercise, we first consider the spread of 
an initial density profile that is concentrated on a ring of 
radius R, 



<r,t = 0) = ^5{r-R) 



(5) 



This example illuminates some pitfalls of the mapping 
convection-diffusion in two dimensions onto a pure dif- 
fusion problem and shows that the size of the ring is 
asymptotically relevant only for sink fiows. 

It is possible to find the spread of the initial ring by 
solving the convection-diffusion equation in two dimen- 
sions subject to the initial condition The solution 
can be written in terms of Besscl functions in the Laplace 
domain or as an infinite series in the time domain. How- 
ever, physical insight about the solution is more easily 
obtained from its asymptotic behavior. To determine 
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the asymptotics, we make use of the mapping ^ and 
express the original problem as a purely diffusive system 
and then use the well-known scaling solution of the latter 
equation to describe the spread from a point mass at the 
origin. 

There is a subtlety in specifying the initial condition. 
The initial two-dimensional ring distribution jSJl is nor- 
malized to unit "mass" 



M{t = 0) = 



z{r, t = 0) 27rr dr = 1 . 



In the effective d-dimensional space, however, the mass 
Mcs is different, 



Mr: 



c{r,t = 0)ndr'^-^dr 



r(i-A) 



(6) 



Here fid = 2tt'^^^ /r{d/2) is the surface area of unit sphere 
in d dimensions. Asymptotically, the concentration pro- 
file should depend on the total effective mass Mcs, while 
the initial size R should become irrelevant. Thus one 
might anticipate that the density approaches 



c(r,t) 



exp 



r 

'4Dt 



(7) 



in the long-time limit. 

Interestingly, Eq. (0 is correct only for sink flows. An 
inconsistency in the ansatz Q) can immediately be seen 
by computing the total mass, M(t) — J c{r,t)2TTr dr, in 
the physical two-dimensional space. Using (U)) we obtain 



M(i) 



1 /^^ZJiV 



(8) 



in the long-time limit. Equation ||SJ) looks reasonable for 
sink flows — the sink causes particles to disappear which 
is reflected by the decay of M(t). On the other hand, 
for source flows (A > 0) equation ||SJ) is invalid — the pre- 
diction that the total mass diverges as is clearly non- 
sensical. 

To resolve this puzzle, it is useful to recall the deriva- 
tion of Eq. (TJ to see why it does not apply to source 
flows. The fundamental solution Q may be found by 
noting that the diffusion equation is invariant under the 
transformation r — > ar, t — > a^t. Hence the scaling vari- 
able r] = /4Dt remains invariant under this scale trans- 
formation, suggesting that the fundamental solution has 
the form c{r,t) = t~^F{r]). Then mass conservation 
/ d'^rc{r,t) — M — const, yields h — d/2. Finally, by 
substituting c(r, t) — t^'^/'^F{rf) into the diffusion equa- 
tion and solving the resulting ordinary differential equa- 
tion we find F = e"'' and thus recover Eq. 

For the radial flow problem, however, mass conserva- 
tion arises in the physical two-dimensional space, rather 
than in the effective d-dimensional space. Therefore, 



b = dpiiys/2 = 1; that is, the asymptotic solution should 
read 



(9) 



SO that the conservation law 



/•OO pOO 

M{t)^ c{r,t)2Trrdr = F{rj)dri = l (10) 
Jo Jo 

indeed holds. By substituting into the governing dif- 
fusion equation (j^J we find 



rjF" + {1~X + Tj)F' + F = 0, 
which is solved to yield 



r(i + A) 



(11) 



(12) 



Thus for source flows, the total mass remains constant, 
M{t) = 1, and the density is 



c(r, t) 



1 1 f Y f_jl\ 

r(l + A) AnDt \4Dt J ^^"^ \~ ADt J 



In this case, the system "forgets" the initial conditions. 

Conversely, our hand-waving argument is correct for 
sink flows as can be verified by explicit solution of the 
full convection-diffusion equation |y| . Therefore the total 
mass decays according to Eq. ((SJ and the density is given 
by Q, or equivalcntly 



c(r, t) 



1 



1 



(—)' 

r(l - A) AirDt \ R"^ J 



cxp 



ADt) 



The system "remembers" initial conditions as it is re- 
flected by the appearance of R in the above asymptotic 
expressions. 



III. DISTANCE TO THE CLOSEST PARTICLE 

As an illustration of the utility of the equivalence be- 
tween convection-diffusion in two dimensions and pure 
diffusion in general dimensions, we examine here the fol- 
lowing question: What is the typical distance from the 
boundary of the absorbing trap to the closest particle? 
Suppose that the trap is a disc of radius R centered at 
the origin, so that particles freely diffuse for r > R and 
get absorbed aX r — R. We assume that initial density is 
uniform, c\t=o = cq. We anticipate that the typical sep- 
aration Tniin between the trap and the closest surviving 
particle exhibits the following behaviors: 

• For sink flows, rmin is finite. 



• For source flows. 



diverges with time. 
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We now confirm these expectations and establish the pre- 
cise asymptotic form of the minimal distance rmin by em- 
ploying a simple quasi-static approach 0, 13| ■ The mini- 
mal distance can also be found by more formal methods 

For sink flows, the concentration approaches a time- 
independent profile Coo(r) that is found by solving the 
Laplace equation subject to the absorbing boundary con- 
dition, Coo\r=R — 0. This solution is 



Af) = Co 



(13) 



We now estimate rmin from the extreme- value criterion 



c(r, t)2'Kr dr — 1 . 



(14) 



which for the density profile (|13|l becomes 

1 = J Co 1 - f — j 2nrdr. 

Computing the integral we find that the dimcnsionlcss 
minimal distance Mmin = f^in/R obeys 



,2+2A 



2^-1 = 



- 1 - 



1 



1 + A 



For source fiows, the minimal distance grows with time 
because particles are driven away from the origin. In us- 
ing the connection to pure diffusion in a space of dimen- 
sion dcff, we note that d^s < 2 and therefore we can set 
the radius of the trap to zero because a diffusing parti- 
cle is recurrent and will hit a point trap with certainty 
if doff < 2. This equivalent problem has the advantage 
of being tractable analytically by elementary means. In- 
deed, we look for a scaling solution of the form 



c(r,t) = coF(?7), V^^, 



that also satisfies 



C\r=0 = 0, c\t^o = Cq. 



(15) 



(16) 



By substituting the scaling ansatz (|15|l into the governing 
diffusion equation (PJ we find 



r,F" + (1 - A + J])F' = 0, 
whose solution is 



dC 



A-l 



r(A) 



(17) 



(18) 



Using Eq. (|14|) together with (|18ll . the minimal value 
Vmin = of tlic scafiug variable obeys 



r(A) 



In the long-time limit, r]„ 



and hence the above 



integral simplifies to 'r\'A+2^)"' ^mln • Thus we arrive at 



r(A 



TTCo 



{ADty 



l/(2A+2) 



(19) 



Let us now illustrate our results by considering some 
specific examples. For sink flow with A = —1/2, the 
effective dimension dcff — 3 and 

r-min = i?+ (Trco)-^/^ 

For source flow with A = 1/2, dcff = 1 and 



9 Dt 



1/6 



For source flow with A = 1, doff = and 

,7r Co y 

Note that for true one-dimensional diffusion, the distance 
from the trap to the closest surviving particle also scales 
as {Dt/ciy/'^ Finally for very strong source 

flows with A ;s> 1 , the effective dimension is d <C and 



In this case, the particles are driven away from the trap in 
an essentially deterministic manner and the above scaling 
dependence of rmin follows directly from r = v ^ Q/f. 

The case of pure diffusion in two dimensions is most 
subtle: the density profile remains time-dependent, but 
we cannot seek a solution in a scaling form because the 
size of the trap cannot be ignored in two dimensions . 
The divergence of the minimal distance with time is also 
particularly delicate in this case, namely 01^3) 



\n{Dt/R^) 



2ttco ln[ln(i:it/i?2)] 



(20) 



IV. 



INTERACTING PARTICLE SYSTEMS 
WITH A LOCALIZED SOURCE 



Consider now the influence of radial flow on an in- 
teracting particle systems. Particularly simple models 
of this type are irreversible annihilation or irreversible 
coalescence 0, 13, 14J. Suppose that the system 
is initially empty and that a localized particle source of 
strength J at the origin is turned on at time t = 0. Par- 
ticles are emitted at the origin and are advected by the 
radial flow fleld. The particles also diffuse and react upon 
colliding with each other (we set the particle radii equal 
to one). These processes may be described by the rate 
equation 



dc / 



-tcfT 



1 dc 
dr 



ln(l/c) 



+ J6{r). (21) 
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Here we employ a phenomenologically-motivated form for 
the reaction term that reflects the true two-dimensional 
spatial nature of the problem. In the homogeneous 
source-free case Eq. (I21f) reduces to c = — Z?c^/ ln(l/c), 
which leads to c{t) ^ [Dt)^^ hi{Dt) asymptotic behavior 
that was rigorously proved to be correct Il5l | . 

We now einploy the same line of reasoning as that given 
in Refs. [l^, [lJ| to determine the particle density. When 
the effective dimension sufficiently large, d > 4, particles 
do not "see" each other; i.e., they interact sufficiently 
weakly that far from the source a non-interacting den- 
sity profile, c(r) ^ ,,-(^-2)^ should arise. For ci < 4, the 
reaction and diffusion terms balance each other, lead- 
ing to c(r) ~ r^^lnr. The borderline case of c? = 4 
needs to be analyzed separately and after some alge- 
bra [see Eqs. (|30|I - H31|I below] it is possible to find a 
nested logarithmic correction to the basic r^^ behavior: 
c(r) ~ r~^[ln(lnr)]^-^. 

Having determined the density profile, we can now 
probe the temporal behavior of the total number of par- 
ticles. Using the fact that the flow field decays as l/r 
and that r = v, the maximum distance traveled by par- 
ticles grows as \/t. Thus the total number of particle 

in the system scales as N{t) ~ drrc{r). Using this 
relation, we arrive at the following conclusions: 



• For sink flows with A < — 1, the effective dimension 
obeys dcs > 4. The density proflle and the total 
number of particles are 



c(r) - r-^''-^) ^ ^2A^ ^(^^ _ ^ 



(22) 



The sink is thus strong enough to keep the total 
number of particles finite. 

• For the case of sink fiow with A = — 1, the effective 
dimension is critical, doff = 4. The density profile 
and the total number of particles are 



c(r) ~ r"2[ln(lnr)]"\ N{t) 



In(lnt)' 



(23) 



The total number of particles therefore exhibits an 
unusually slow growth in this marginal regime. 

• For sink flows with — 1 < A < and also for source 
flows, the effective dimension obeys dcs < 4. The 
density proflle and the total number of particles 
become universal 



c(r) -r"2 In r, 7V(t) - (Ini)^ 



(24) 



We conclude that for the case A > — 1 the total number 
of particles grows extremely slowly with time, while for 
A < — 1 there are just a few particles in the system. 

The results for source flow are seemingly peculiar; here 
one might naively expect that the total density grows 
linearly with time because particles are driven away from 
each other by the flow and thus would not interact. Let us 
therefore re-derive some of our results without exploiting 



the mapping between radial flow in two dimensions to 
pure diffusion in variable spatial dimension. Consider, 
for simplicity, the case of strong source flow, i.e., the 
Peclet number is large. We may then drop the diffusion 
terms and keep just the dominant convection term. Using 
this simplification and assuming that the system is in the 
steady-state regime, the governing equation H21|) becomes 



2A dc 



r dr ln(l/c) 
Integrating this equation we obtain 

ln(l/c) ^ ,^r^ 
^~-^+4A' 



(25) 



(26) 



where we use the shorthand notation A = CQ^ln(l/co). 
We estimate the density cq near the origin from the con- 
servation law 27rrco = J to give cq = J/Q and there- 
fore 



From Eq. (|26|l . we finally obtain the steady state density 
profile 



In 



c(r) 



A 



(27) 



4A 



Sufficiently far from the origin, r :$> VXA, Eq. l(?7| sim- 
plifies to c(r) ~ r^^ Inr in agreement with our previous 
flndings in (|24|l . 

To determine the total number of particles we note 
that t he stea dy state solution formally holds as long as 
r <^ \fQtpK. Using the fact that the leading edge of the 
particles advance as r = ^ Qt/ir, the total number of 
particles density then scales as 



N{t) 



c(r) 2-Kr dr. 



From (EZI) we then obtain 



N{t) - 27rA {\n^{Dt + A) - In^ A] 



(28) 



In the long time limit, Eq. (|28|l qualitatively agrees with 
our previous result, N{t) ^ (Ini)'^. Note, however, that 
logarithmic behavior arises only after a short-time linear 
regime. For Dt < A, we get N{t) - 47rA Dt, which 
simplifies to N{t) ^ Jt. Thus indeed the total density 
initially grows linearly in time in accordance with intu- 
ition. The crossover between these two regimes occurs at 
tc = A/D - (A/J)ln(Q/J). 

The previous analysis equally applies to annihilation 
and coalescence. More generally, we may also consider 
irreversible mass-conserving aggregation where the reac- 
tion of a cluster of mass i with a cluster of mass j leads 
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to a cluster of mass k = i + j. For aggregation, the fun- 
damental quantities are the densities of various particle 
species. Let us set the mass of particles that are emit- 
ted at the origin to unity. Denote the density of these 
monomer particles as ci. Generally let Cfc be the density 
of particles composed of k monomers. We shall seek only 
the stationary densities Cfc(r). These densities satisfy a 
modified Smoluchowski equation 



J'cS 



1 d 



1 



E 



2cfeC 



(29) 



Here we again employed a phenomenologically-motivated 
form for the reaction term (see e.g., |l6j |') and we addition- 
ally assumed that the coalescence rate is independent of 
the masses of the two reactants. (The latter assumption 
was originally made by Smoluchowski; it greatly simpli- 
fies the analysis of the infinite system of rate equations 
[TtI . Furthermore, for diffusion-controlled reactions in 
two dimensions, the reaction rate depends on the reac- 
tant radii in a weak logarithmic fashion 1, 6J.) 

The asymptotic behavior is especially interesting in the 
marginal case of A = — 1 when the effective dimension is 
critical, deff = 4. Then the total density satisfies 



d^c _^ 3 dc 
dr"^ 



r dr ln(l/c) ' 
whose leading asymptotic behavior is 

4 1 



ln(ln r) 

The density of monomers then satisfies 



(30) 



(31) 



dr^ 



3 dci 
r dr 



2c 



Hi/c] 



Cl 



which is a long-time limit of a system of equations orig- 
inally solved by Smoluchowski, see |l7j . with p playing 
the role of time. The solution admits the scaling form, 

Fk = p^^e^'^/P, and therefore 



Cfc 



16 



1 



[In(lnr)]^ 



exp 



4fc 



In(lnr) 



(34) 



Using this result we can verify that the total particle 
density c(r) =J2k>i Cfe(^) is consistent with our previous 
prediction (|31|l : we can also compute the mass density 
m(r) = X]fe>i '^Cfe(r) = r^^ and estimate the total mass 
to be 



M{t) 



m{r) 2'Krdr = tt In t. 



(35) 



When A > —1, that is, the effective dimension is 
smaller than critical, dcff < 4, the particles do not un- 
dergo a sufficient number of collisions to develop a scal- 
ing form like H34|l . Nevertheless, at least the monomer 
density always exhibits interesting dynamical behavior. 
Indeed, using H29() and c = 8(1 + A)r^^ Inr, which is the 
precise form of (|24|l that includes the correct amplitude, 
we find that the monomer density satisfies 



d^ci 1 - 2A dci 8(1 + A) 



dr'^ 

Solving (EH) gives 



dr 



■ Cl 



-A+ v/A2 + 8(l + A), 



(36) 



(37) 



so that the monomer density has a pure algebraic tail 
with a A-dependent exponent. 



V. CONCLUDING REMARKS 



r2(lnr) In (Inr) 



Cl, 



whose solution is 



Cl ^ 



[ln(lnr)]2 



(32) 



Solving for the first few fc-mer densities one by one, we 
observe that they all have the form 



1 



Cfc 



Fk{p), P = 7 In(lnr). 



(33) 



This approach also gives Fk ~ p^^ for small k, but it is 
better to derive this result in general rather than guessing 
this behavior based on the small-/c behavior only. Thus 
substituting the ansatz (|33|) into H29() and keeping only 
the dominant terms we obtain 



m 

dp 



i+]=k ^ 



Allowing the spatial dimension d to be a free parame- 
ter and then developing theoretical approaches based on 
this parameter (an expansion about a critical dimension, 
dimensional regularization, etc.) has proven extremely 
fruitful in field theory and statistical physics 0, Il9| . 
Thus in many treatments of critical phenomena, non- 
integer dimensions arise naturally, but this construction 
is generally used as an intermediate step toward the ulti- 
mate goal of obtaining information in physically relevant 
spatial dimensions, such as d — 1,2,3. However, cer- 
tain complex problems can be mapped onto simpler ones 
that are defined in spaces of a non-integer spatial di- 
mension. We demonstrated that such a connection arises 
for a class of two-dimensional convection-diffusion prob- 
lems in which the flow field is radial and proportional 
to 1/r. It would be exciting to realize these convection- 
diffusion fiows experimentally and hence probe dynami- 
cal processes in spaces whose dimension is not necessarily 
an integer. 
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